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Microprocessor  Requirements  for  Implementing 
Modern  Control  Logic 


SUMMARY 


Analytical  procedures  for  establishing  microprocessor  requirements  for 
multivariable  feedback  control  of  linear  stochastic  dynamic  systems  were 
developed  and  evaluated.  Key  Issues  In  microprocessor  Implementation  of  modern 
estimation  and  control  logic  Include  (1)  accuracy,  (2)  computational  capability 
including  arithmetic  as  well  as  Interface  speed,  and  (3)  memory  requirements. 
Compatibility  of  current  microprocessor  technology  with  the  established  implemen- 
tation requirements  was  assessed. 

The  analytical  method  to  establish  microprocessor  accuracy  requirements  for 
closed-loop  digital  control  of  linear  stochastic  systems  is  based  upon  defining 
a quadratic  performance  index  that  provides  a quantitative  measure  of  performance 
degradation.  This  method  Is  employed  to  determine  closed-loop  system  time 
response  sensitivity  to  (1)  microprocessor  word  length  and  (2)  system  structure 
(standard,  Jordan  canonical  or  companion  forms)  within  the  control  logic.  The 
technique  for  determining  computational  capability  requirements  Include  proce- 
dures for  (1)  defining  the  maximum  controller  sample  time  as  a function  of 
system  dynamics  and  (2)  calculating  controller  arithmetic  and  Interface  computa- 
tion time  per  sampling  Interval.  Memory  requirements  are  delineated  as  a 
function  of  system  model  and  microprocessor. 

The  developed  procedures  were  evaluated  and  illustrated  by  application  to 
(1)  a second-order  system  and  (2)  a linearized  fifth-order  F100  turbofan  engine 
model.  Results  were  verified  using  a digital  simulation  of  a continuous  system/ 
microprocessor  controller.  Microprocessors  considered  were  the  Intel  8080  micro- 
processor (an  8 bit  microprocessor)  and  the  Digital  Equipment  Corporation  (DEC) 
LSI  11/2  microprocessor  (a  16  bit  microprocessor) . 
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RESULTS  AND  CONCLUSIONS 


1.  Digital  electronic  requirements  for  Implementing  optimal  stochastic 
feedback  controls  designed  using  linear  quadratic  Gaussian  (LQG)  theory  were 
delineated.  These  requirements  Include  accuracy,  computational  capability,  and 
memory  requirements.  The  effect  on  these  requirements  of  various  system  represen- 
tations (standard,  Jordan  canonical  and  companion  forms)  within  the  control  logic 
was  investigated. 

2.  Accuracy  requirements  depend  upon  closed-loop  system  dynamics.  To 
establish  accuracy  requirements  a quadratic  performance  index  was  defined  to 
provide  a direct  indication  of  transient  performance  degradation  associated  with 
finite  microprocessor  word  lengths.  Results  show  that  performance  degradation 
as  a function  of  word  length  depends  upon  model  structure  within  the  controller. 
Therefore,  the  performance  index  approach  is  used  to  determine  the  model  structure 
to  minimize  performance  degradation  as  well  as  to  determine  microprocessor  word 
length. 


3.  Computational  requirements  depend  upon  (1)  the  maximum  sample  time  and 
(2)  controller  computation  time  per  sampling  interval.  The  computation  time  per 
sampling  interval  must  be  less  than  the  maximum  sample  time. 

(a)  The  maximum  sample  time  is  a function  of  closed-loop  system  dynamics. 
Plotting  the  z-plane  root  locus  of  the  discrete  controller  as  a function 
of  sample  time  was  shown  in  this  study  to  be  an  effective  way  to 
determine  the  effect  on  system  performance  of  different  sample  times. 

By  analyzing  the  z-plane  root-locus  plot  the  maximum  sample  time  is 
established. 

(b)  The  computation  time  per  sampling  interval  includes  arithmetic  and  I/O 
computation  times.  The  arithmetic  computation  time  is  a function  of 
the  number  of  system  states,  inputs,  and  measurements  as  well  as  micro- 
processor speed.  The  I/O  computation  time  depends  upon  the  number 

of  system  inputs  and  measurements  as  well  as  interface  speed.  Computa- 
tion times  for  the  Intel  8080  microprocessor  (an  8 bit  microprocessor), 
the  DEC  LSI  11/2  microprocessor  (a  16  bit  microprocessor) , and  the 
DEC  LSI  11/2  microprocessor  with  hardware  multiply  option  were  tabu- 
lated for  the  generic  linear  system  as  a function  of  system  order, 
number  of  inputs  and  outputs.  These  tables  indicate  that  the  DEC  LSI 
11/2  microprocessor  with  hardware  multiply  is  significantly  faster 
than  the  Intel  8060  and  DEC  LSI  11/2  microprocessors.  For  implementing 
LQG  control  logic  the  Intel  8080  microprocessor  is  faster  than  the  DEC 
LSI  11/2  microprocessor  because  the  Intel  8080  software  multiply 
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subroutine  (8  bits)  Is  significantly  faster  than  the  DEC  LSI  11/2 
software  multiply  subroutine  (16  bits).  In  addition,  the  Jordan 
canonical  structure  requires  the  minimum  arithmetic  computation 
time;  whereas,  the  standard  structure  requires  the  maximum  arithmetic 
computation  time. 

4.  Memory  requirements  depend  upon  the  system  model  and  microprocessor  code. 
Memory  requirements  for  the  Jordan  canonical  and  companion  structures  are  the 
same.  The  standard  structure  requires  the  maximum  memory.  Memory  requirements 
for  the  Intel  8080  and  DEC  LSI  11/2  microprocessors  were  tabulated.  These  tables 
show  that  memory  requirements  for  LQG  control  logic  are  not  severe.  For  example, 
a controller  for  a system  with  10  states,  10  inputs  and  10  outputs  requires 

less  than  IK  of  RAM  and  less  than  IK  of  PROM. 

5.  The  developed  procedures  were  validated  by  application  to  two  systems. 
The  first  system  was  a single- input,  single-measurement  second-order  linear 
system.  The  second  plant  was  an  F100  gas  turbine  engine  model  linearized  about 
sea-level  static  military  operation.  The  engine  model  has  five  states,  four 
inputs,  and  five  measurements.  The  procedures  were  verified  using  a closed- 
loop  simulation  of  the  continuous  plant/discrete  controller. 

(a)  An  Intel  8080  microprocessor  can  be  employed  to  Implement  the  second- 
order  controller.  The  requirements  for  the  second-order  controller 
are  (1)  8-bit  word  length  accuracy,  (2)  minimum  and  maximum  sample 
times  of  2.73  and  700  msec,  respectively,  and  (3)  347  and  21  words  of 
PROM  and  RAM,  respectively. 

(b)  A DEC  LSI  11/2  microprocessor  with  hardware  multiply  can  be  used  to 
implement  the  LQG  F100  engine  controller.  The  requirements  for  the 
F100  engine  controller  are  (1)  12-bit  word  length  accuracy,  (2)  minimum 
and  maximum  sample  times  of  9.68  and  25.0  msec,  respectively,  and 

(3)  less  than  200  words  of  RAM  and  less  than  200  words  of  PROM. 
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INTRODUCTION 


Over  the  past  several  years  use  of  modern  control  methodology  — In  particular, 
linear  quadratic  Gaussian  (LQG)  theory  — has  gained  Increased  recognition  as  an 
effective  design  tool  for  control  of  nonlinear  multivariable  stochastic  systems 
(Refs.  1-8).  The  referenced  studies  have  been  conducted  under  a combination  of 
AFOSR,  Office  of  Naval  Research  (ONR) , Air  Force  Aero  Propulsion  Laboratory  (AFAPL) , 
NASA-Lewis  and  Pratt  & Whitney  Aircraft  (P&WA)  support.  In  these  as  well  as  many 
other  aerospace  applications  the  primary  impetus  for  application  of  modern  LQG 
control  concepts  is  Improved  system  performance  combined  with  the  advent  of  digital 
electronic  control  Implementation.  Digital  electronics  provide  the  means  by 
which  complex  controllers  associated  with  LQG  theory  can  be  Implemented.  The 
current  trend  both  within  as  well  as  outside  the  aerospace  controls  community 
toward  increased  use  of  digital  electronics  — in  particular,  microprocessors  — 
will  lead  to  increased  use  of  modern  control  logic  including  system  identification, 
modeling,  estimation,  and  multivariable  control  methodologies  (Ref.  9).  In  addi- 
tion, use  of  microcomputer  controllers  will  lead,  in  many  instances,  to  reduced 
control  cost  (Refs.  10  and  11),  lighter  and  smaller  controls  (Ref.  12),  lower 
power  requirements  and  integrated  circuit  reliability  (Ref.  13).  Recent  studies 
(Ref.  14)  have  demonstrated  that  existing  microprocessors  can  be  used  to  imple- 
ment algorithms  for  parameter  identification  of  relatively  simple,  low-order 
dynamic  systems. 

However,  prior  to  widespread  use  of  microprocessors  for  modern  control  logic 
Implementation,  key  issues  associated  with  microprocessor  Implementation  of  LQG 
control  and  estimation  concepts  must  be  addressed  and  resolved.  These  issues 
Include  (1)  accuracy,  (2)  computational  capability  including  arithmetic  as  well 
as  interface  speed,  and  (3)  memory  requirements  (Ref.  15).  These  requirements 
depend  upon  system  dynamics  as  well  as  upon  the  particular  control  algorithm 
employed.  Definin''  these  requirements  will  establish  criteria  for  selecting 
the  appropriate  computer  system  for  control  implementation. 

Consequently,  this  study  was  directed  toward  establishing  microprocessor 
requirements  for  implementing  modern  control  logic.  The  objective  of  this 
investigation  was  to  define  microprocessor  requirements  of  LQG  control  logic 
for  linear  systems.  An  analytical  procedure  to  determine  accuracy  requirements 
for  LQG  control  logic  was  developed.  Computational  capability,  interface  and 
memory  requirements  were  established  for  generic  representations  of  system 
dynamics.  The  effects  on  these  requirements  of  various  system  representations 
(standard,  Jordan  canonical,  and  companion  forms)  within  the  LQG  control  logic 
were  investigated.  To  evaluate  and  illustrate  the  developed  methodology,  micro- 
processor LQG  control  requirements  were  delineated  for  (1)  a second-order  model 
and  (2)  a fifth-order  F100  turbofan  engine  model. 
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Because  the  linear  system  model  employed  In  this  study  Is  general  in  nature, 
results  of  this  study  are  not  unique  to  any  one  particular  system.  Therefore, 
the  results  obtained  in  this  study  have  applicability  to  a broad  range  of  control 
problems  with  which  the  Air  Force  and  the  other  military  services  are  concerned. 
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CONTROL  OF  LINEAR  STOCHASTIC  SYSTEMS 


In  this  section  the  linear  stochastic  control  problem  is  presented.  The 
linear  stochastic  system  model  is  defined  first.  This  model  is  general  in 
nature  and  is  therefore  applicable  to  a broad  range  of  estimation  and  control 
problems.  Linear  quadratic  Gaussian  (LQG)  control  design  is  then  outlined.  In 
the  final  part  of  this  section  digital  implementation  of  the  multivariable  control 
and  estimation  logic  is  discussed. 


System  Description 

The  system  model  is  shown  in  Fig.  1 with  provision  for  estimation  and  control 
algorithms  included.  The  open-loop  system  consists  of  a linear  plant,  control 
actuators,  and  sensois.  Open-loop  linear  system  dynamics  are  described  by  the 
differential  and  algebraic  equations 

X(t)  = A X(t)+  Bu(t)  + £(t) 
y(t)  - c x (t)  + D u(t) 
z(t)  = E x(t)+rj(t) 

where  the  vectors  x,  u,  y,  and  z represent  the  n states,  m inputs,  p outputs,  and 
£ measurements,  respectively.  Note  that  the  vector  x includes  plant,  actuator, 
and  sensor  states.  The  random  process  vectors  £ and  q represent  white  zero-mean 
Gaussian  n-dimensional  process  and  ^-dimensional  measurement  noise,  respectively. 
The  constant  A (n  x n) , B (n  x m)  , C (p  x n) , D (p  x m)  and  E ( £ x n)  matrices 
define  linear  time-invariant  system  dynamics.  The  dot  notation  denotes  differen- 
tiation with  respect  to  time.  The  Initial  state  vector  is  assumed  to  be  a Gaussian 
random  variable  with  known  mean.  The  random  vectors  x(0) , £(t),  and  n ( t)  are 
assumed  independent  with  known  covariances.  The  statistical  properties  of  these 
random  vectors  define  the  system  inaccuracies.  The  process  noise  (£(t)) 
models  actuator  uncertainties,  plant  disturbances,  and  system-to-system  parameter 
variations.  Sensor  noise  (n(t))  models  measurement  inaccuracies. 

Linear  Quadratic  Gaussian  Control 

Linear  quadratic  Gaussian  theory  provides  the  analytical  tools  to  design  an 
optimal  feedback  control  for  the  stochastic  system  described  by  Eq.  (1).  Solution 
of  the  overall  LQG  problem  separates  into  the  solution  of  (1)  linear  quadratic 
regulator  (LQR)  and  (2)  linear  Gaussian  estimation  problems.  The  key  theorem 
that  demonstrates  this  property  is  often  referred  to  as  the  separation  theorem 
(Ref.  16).  The  stochastic  control  design  procedures  then  involve  (1)  determin- 
istic multivariable  feedback  control  using  LQR  theory, (2)  stochastic  filter  design 
using  Kalman  estimation  logic,  and  (3)  closed-loop  regulation  based  on  feedback 
of  estimated  system  variables  through  the  deterministic  control  logic. 
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The  solution  to  the  stochastic  linear  optimal  control  problem  is  well-known 
and  the  details  can  be  found  in  Ref.  16.  Under  appropriate  controllability  and 
observability  conditions  the  optimal  input  u*  exists.  The  filter  and  control 
dynamics  are  described  by  *(t)=  F*(t)  + H(Z(t) -E$(t)) 

U*(t)  = Gft(t)  (2) 

F= A+BG 

where  the  notation  (A)  denotes  the  estimate  of  the  variable  in  parentheses, 

G (m  x n)  represents  the  deterministic  feedback  control  gain  matrix  and  H (n  x t ) 
represents  the  steady-state  Kalman  filter  gain  matrix. 

The  deterministic  control  and  Kalman  filter  gain  matrices  (G  and  H,  respec- 
tively) depend  upon  (1)  the  dynamics  of  the  system,  (2)  the  levels  of  uncertainty 
in  the  system,  and  (3)  performance  criteria  that  specify  satisfactory  time  evo- 
lution of  the  system  inputs  and  outputs.  However,  since  the  problem  assumptions 
lead  to  constant  gain  matrices,  these  gains  can  be  computed  off-line.  The  result- 
ing stochastic  feedback  control  structure  to  be  implemented  using  microprocessors 
is  illustrated  in  Fig.  2. 


Digital  Control  Implementation 


The  structure  of  the  stochastic  feedback  control  to  be  implemented  with  micro- 
processor is  shown  in  Fig.  2.  Figure  2 indicates  that  a system  model  must  be 
employed  within  the  controller.  Many  mathematically  equivalent  model  representa- 
tions of  the  linear  system  dynamics  exist  (Ref.  17).  The  model  representation 
employed  will  affect  the  digital  electronic  requirements.  Generally,  for  physical 
systems  the  matrix  F will  be  standard  form.  For  the  standard  form  all  elements 
of  the  matrix  F are  assumed  to  be  nonzero  and  F is  represented  by 


f||  *12  * * f|n 

^ 21  ^ 22  ’ * ^2n 


(3) 


F = 


fnl  W 


'nn 


To  reduce  digital  electronic  requirements  for  implementing  the  LQG  control 
logic  (Eq.  (2)),  state  coordinates  employed  within  the  controller  may  be  different 
from  the  physical  system  state  coordinates.  Several  forms  obtained  by  state 
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transformation  (i.e.,  by  a change  of  coordinates)  will  reduce  the  required  number 
of  arithmetic  operations.  These  forms  include  the  Jordan  canonical  form  and 
the  companion  form. 

To  change  state  coordinates  the  estimated  state  vector  x is  transformed 
through  the  equation 

A _-i  A 

W = T 1 X 

(4) 

where  T is  an  n x n nonsingular  constant  matrix.  Substituting  Eq.  (4)  into  the 
state  equations  for  x (Eq.  (2))  leads  to  transformed  state  equations  given  by 

fc(t)  = T''FT$(t)+  T'1  H[z(t)-  ETd(t)]  ( 

u(t)  = 6T0(t). 

Equations  (2)  and  (5)  represent  the  same  system  in  different  coordinates.  Note, 
also,  that  Eq.  (5)  is  of  the  same  form  as  Eq.  (2).  That  is, 

4(t)  = FTw(t)+HT(z(t)-ETwlt))  (6) 

u(t)=  Gtw  (t) 

where  FT  - T-1(F)T,  Up  - T_1H,  ET  - ET  and  GT  ■ GT.  However,  the  matrix  FT  depends 
upon  the  selected  transformation  matrix  T.  Note  that  T-I  results  in  the  standard 
form,  i.e.,  FT“F. 

To  transform  Eq.  (2)  to  the  Jordan  canonical  form,  the  eigenvalues  of  the 
matrix  F must  be  employed.  For  convenience  it  is  assumed  here  that  the  matrix 
F possesses  distinct  eigenvalues.  This  assumption  is  quite  reasonable  for  most 
physical  systems.  Associated  with  each  eigenvalue  X^  there  exists  a nonzero 
eigenvector  (n  x 1)  defined  by 

F u j = A j u)  j , /7\ 


When  the  eigenvalues  are  distinct  the  eigenvectors  form  a linearly  independent 
set.  The  matrix  T to  transform  Eq.  (2)  to  the  Jordan  canonical  form  is  the  non- 
singular modal  matrix  of  F (i.e.,  the  columns  of  T are  the  eigenvectors  u^) • The 
matrix  T~1  FT  is  a diagonal  matrix  with  diagonal  elements  equal  to  the  eigenvalues 
of  F,  i.e. , 

Tx,  o • • o ! 


T FT* 
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where  T Is  the  modal  matrix  of  F. 

To  transform  the  matrix  F to  the  companion  form,  the  matrix  T (which  Is  not 
unique  (Ref.  16))  may  be  selected  as 


where 


h(l  xn)  = (1,0, ,0  ]. 


The  matrix  T_1FT  is  then  the  companion  form  given  by 


T“'FT 


Qq  Qi  02 
a"  On  On 


-2nd 


1 1 

I s 


where  a^,  1-0,  1..,  n are  the  coefficients  of  the  characteristic  polynomial  of  F. 
That  Is, 


an  Xn  + a„_,  + o0=|f-\i|. 


Inspection  of  the  standard,  Jordan  canonical,  and  companion  forms  (Gqs.  (3), 
(8)  and  (11),  respectively)  and  the  filter  equation  (Eq.  (5))  Indicates  that  the 
required  number  of  arithmetic  operations  (additions  and  multiplications)  is  a 
function  of  model  representation.  The  model  representation  selected  on  the  basis 
of  minimum  number  of  arithmetic  operations  would  be  the  Jordan  canonical  form. 
However,  the  Ref.  18  study  shows  that  accuracy  of  digital  filters  varies  with 
filter  structure  even  though  the  overall  transfer  functions  of  the  digital  filters 
are  Identical.  This  variance  in  accuracy  is  due  to  the  finite  word  lengths 
associated  with  digital  electronic  implementation.  Therefore,  this  study 
investigated  the  impact  of  filter  structure  on  system  performance.  The  approach 
adopted  for  determining  the  appropriate  structure  is  discussed  in  the  next 
section. 
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To  code  the  controller  on  a digital  computer  the  filter  equations  for  state 
estimation  (Eq.  (5))  are  generally  represented  by  (1)  state  prediction  equations 
and  (2)  state  update  equations.  Filter  equations  — obtained  by  expanding  closed- 
loop  filter  dynamics  in  a Taylor  series  — which  predict  state  variables  at  time 
(t  + At),  given  measurements  to  time  t,  are  described  by 

$(t+At/t)  = (i  + *)$(t/t) 


4>  = T*1  FT  At  + 


(T'1  FT)2  At*  . (T*1  FT)*  At® 


When  the  sampling  interval  has  been  chosen,  $ may  be  computed  off-line  to  the 
desired  accuracy.  The  number  of  terms  retained  in  the  series  expansion  for  $ 
depends  upon  the  microprocessor  word  length  to  be  used  in  implementing  the  con- 
troller. That  is,  the  number  of  terms  retained  is  selected  so  that  the  sum  of 
terms  beyond  the  last  term  retained  is  less  than  the  accuracy  of  the  microprocessor. 
The  notation  f(t  + At/t)  represents  the  value  of  the  function  f at  tlme(t  + At) 
given  measurements  to  time  t.  The  filter  update  and  deterministic  control  equa- 
tions are  given  by 

w(t  + At/t+At)  *fc(t+At/t)  + r'HAtfztt+AtJ-ET&tt  + At/t)! 

* 1 (14) 

u (t  + At/t  + At)  = GT  G(t  + At/t  + At)  . 


Equation  (13)  indicates  that  state  prediction  computational  requirements  depend 
upon  the  system  model  employed  in  the  filter  and  the  state  order;  whereas,  Eq. 
(14)  shows  tjhat  state  update  and  input  computational  requirements  depend  on  state 
and  measurement  orders  and  input  and  state  orders,  respectively.  In  addition, 
Eqs.  (13)  add  (14)  show  that  the  arithmetic  operations  for  solving  the  filter 
and  deterministic  control  equations  are  matrix  multiplications  and  additions. 


For  digjltal  implementation  system  dynamics  are  normalized  so  that  the  control 
law  (Eqs.  (1|3)  and  (14))  may  be  coded  using  fixed-point  rather  than  floating- 
point arithmetic.  Fixed-point  arithmetic  is  employed  for  digital  implementation 
of  the  control  logic  to  increase  computational  speed.  Generally,  system  dynamics 
are  normalized  so  that  the  digital  numbers  in  the  control  computations  range 
between  -1.0  land  +1.0.  Then  the  most  significant  bit  (bit  0)  in  the  computer 
word  is  the  sign  bit  and  the  least  significant  bit  (bit  b-1;  where  b denotes  the 
number  of  bit*  in  the  microprocessor  word)  represents  2(_b+1).  For  notational 
convenience,  (1 . e . , to  avoid  introducing  new  symbols  for  normalized  system  vari- 
ables) Eqs.  (l)  through  (14)  will  be  used  to  represent  normalized  system  dynamics 
in  the  remainder  of  this  report. 


Accuracy  of  the  microprocessor  estimation  and  control  logic  will  depend  upon 
(1)  word  length  of  the  microprocessor  and  (2)  word  length  of  the  input/output 
(I/O)  devices.  Due  to  finite  microprocessor  and  I/O  word  lengths  (where  I/O 
devices  include  direct  digital  output  effectors  such  as  transducers  with  frequency 
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output  and  stepper  motors  as  veil  as  A/D  and  D/A  converters)  quantization  errors 
occur  In  the  digital  representation  of  a number . For  the  scaling  described  pre- 
viously, the  absolute  value  of  the  maximum  quantisation  error  e la  given  by 


R7 9-944258-2 


MICROPROCESSOR  REQUIREMENTS 


Key  issues  associated  with  microprocessor  implementation  of  linear  quadratic 
Gaussian  (LQG)  control  logic  are  addressed  in  this  section.  An  analytical  pro- 
cedure to  establish  microprocessor  accuracy  requirements  is  described  first. 
Computational  requirements  including  arithmetic  as  well  as  Interface  speed  are 
then  set  forth.  Memory  requirements  as  a function  of  system  order  and  system 
model  within  the  filter  are  then  presented. 


Accuracy  Requirements 

Controller  accuracy  required  to  obtain  satisfactory  system  performance  will 
depend  upon  (1)  the  structure  of  the  system  model  employed  within  the  Kalman  filter 
and  (2)  system  sensitivity  to  controller  errors.  To  determine  system  time  response 
sensitivity  to  microprocessor  word  length  (that  is,  to  determine  microprocessor 
accuracy  requirements)  a performance  index 


/CD 

[(SYSTEM  RESPONSE) ‘—(SYSTEM  RESPONSE)1]3  dt 


(16) 


— where  ( )*  represents  the  system  response  with  the  controller  coded  on  a 36-bit 

floating-point  computer  (a  very  accurate  controller)  and  ( )+  represents  the 

system  response  with  the  controller  coded  using  b bits  (with  b - 36)  and  fixed- 
point  arithmetic — was  defined  to  provide  a quantitative  measure  of  performance 
degradation.  Its  computed  value  provides  a direct  Indication  of  the  transient 
performance  degradation  associated  with  microprocessor  quantization  errors. 

The  Integrand  of  J was  selected  as  quadratic  in  system  input  and  output 
variables  so  that  the  value  of  J can  be  analytically  computed  from  the  known  pro- 
blem matrices.  That  is,  the  performance  index 


o(y*-yt)  ♦(u’-uty  R(u*-u^)Jdt 


(17) 


where  y*  m output  response  vector  with  36-bit  controller 

y*  • output  response  vector  with  b-bit  controller 
u*  » control  vector  with  36-bit  controller 
u+  - control  vector  with  b-bit  controller 
Q,R  - weighting  matrices 

was  defined.  The  performance  index  J represents  performance  degradation  due  to 
finite  word  lengths  less  than  the  accurate  36-blt  word  length.  As  the  number 
of  bits  in  the  computer  word  approaches  36,  J approaches  zero. 
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From  Eqs.  (1)  and  (5)  the  closed-loop  system  dynamics  (open-loop  system  and 
controller)  are  described  by 


X » F X 

^ a# 

y * C x 
u=S  x 


7(0). 


(18) 


X s 


‘ X " 

~ a“ 

. W . 


where 


F = 


h)ce|  (t'1  FT)c  -(t'1  H)c  (et)c J 

c = [c  j o(gt)c ] 

G = [o  j (GT)C]  • 


I 

I 

■+■ 

I 


b(gt)# 


(19) 


The  notation  ( )c  Indicates  that  the  matrix  is  Implemented  in  the  digital  con- 
troller. The  accuracy  of  these  matrices  in  the  controller  is  limited  by  the  word 
length  of  the  digital  controller.  Matrices  without  the  subscript  c represent  the 
open-loop  system  dynamics.  The  solution  to  Eq.  (18)  is  given  by 

x s ePT  x0 


m £?*  ^ 

u*Ge  x- 


(20) 


Substituting  Eq.  (20)  for  the  36-bit  representation  (y*,  u*)  and  for  the 
b-bit  representation  (y+,  u+)  into  Eq.  (17)  and  algebraically  manipulating  the 
resulting  equation  for  J as  well  as  taking  the  expected  value  of  J over  xQ  leads 
to 


where 


j«  tr  P,  ♦ 2trP2  ♦tr  Ps 

?*'p«»ptF,.(c*'0ct  .*•'**') 

?*'p,.p,  ♦3t'p8t) 


(21) 


(22) 


and  "tr"  indicates  the  trace  of  a square  matrix.  Equations  (21)  and  (22)  indicate 
that  the  performance  index  J is  a function  of  (1)  the  closed-loop  system  dynamics 
(F,  C,  6)  and  (2)  the  microprocessor  word  length  (number  of  bits) . For  a given 
closed-loop  system  the  performance  degradation  J averaged  over  all  possible  initial 
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conditions  can  be  computed  by  solving  the  bilinear  equations  (Eq.  (22))  for 
P.  , P2,  and  P^.  This  analytical  procedure  eliminates  the  need  to  run  system 
time  responses  and  perform  numerical  integrations  over  many  initial  conditions 
to  determine  microprocessor  accuracy  requirements  for  implementing  the  digital 
controller.  In  addition,  this  procedure  for  establishing  microprocessor  accuracy 
requirements  is  not  limited  to  LQG  controllers.  The  procedure  is  applicable  to 
any  linear  closed-loop  system  expressed  in  Eq.  (18)  form. 

In  addition  to  using  Eq.  (22)  to  establish  the  required  microprocessor 
word  length,  Eq.  (22)  can  also  be  employed  to  determine  the  structure  (i.e., 

F — see  Eq.  (19))  to  minimize  performance  degradation.  That  is,  time  response 
y*»  u*  associated  with  the  high  accuracy  36-bit  floating-point  controller  should 
have  negligible  dependence  upon  system  structure  since  the  system  structures  F 
are  mathematically  equivalent.  On  the  other  hand,  time  response  y*,  u^  will 
depend  upon  the  particular  structure  selected  because  of  fixed-point  arithmetic 
operations  and  small  word  length  (e.g.,  b = 4 or  8) . Results  obtained  will 
assist  in  establishing  the  optimum  structure  for  microprocessor  implementation 
of  LQG  control  logic. 


Computational  Requirements 

In  this  section  computational  requirements  are  discussed.  Computational 
requirements  depend  upon  (1)  the  maximum  sample  time  (i.e.,  minimum  sampling 
rate)  that  maintains  desired  system  performance  and  system  stability  as  well  as 
(2)  arithmetic  and  I/O  computation  time  per  sampling  interval.  A procedure  for 
determining  the  maximum  sample  time  (i.e.,  maximum  sampling  interval)  is  described 
first.  The  maximum  sample  time  depends  upon  system  dynamics.  Computation  time 
per  sampling  interval  is  then  presented.  The  computation  time  per  sampling  inter- 
val depends  upon  (1)  the  number  of  arithmetic  and  I/O  operations  per  sampling 
interval,  and  (2)  the  speed  of  the  microprocessor  and  Interface.  Computation 
time  per  sampling  interval  establishes  the  minimum  sample  time.  If  required  com- 
putation time  per  sampling  interval  is  less  than  the  maximum  sample  time  that 
maintains  desired  system  response,  then  the  microprocessor  and  Interface  system 
may  be  used  to  implement  the  LQG  controller. 


Sample  Time 

To  analyze  linear,  time-invariant,  discrete-time  systems  (Eqs.  (13)  and  (14)), 
the  z-transform  is  employed.  From  Eqs.  (13)  and  (14),  discrete  controller  dynsmlcs 
are  given  by 
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where 


d(k  + i)s^0w(k)  + H0  Z (k  + l) 
u(k*i)*  G0  w(k  + i) 


<^0=(l-THHETAt)  (l  + ^>) 


4>  = T*'  FTAt  + 
H0  = T*'HAt 
G0  = GT 


(T*1  FT)2  At 2 + (T~'FT)3At3 


2! 


3! 


(23) 


(24) 


and  k denotes  the  kth  sample  time.  Routh's  criterion,  the  root-locus  method, 
or  frequency  response  methods  — procedures  used  to  investigate  transient  re- 
sponse and  stability  of  continuous-time  systems  — may  be  extended  to  investi- 
gate transient  response  and  stability  of  discrete-time  systems  (Ref.  19). 


In  this  study,  the  root-locus  method  was  employed  to  relate  sample  time 
(At)  and  system  performance.  This  method  has  the  advantage  that  the  z-plane 
poles  are  available  on  the  root-locus  plot  once  the  controller  has  been  selected 
and  the  root  locus  plotted  as  a function  of  sample  time.  The  effect  on  system 
performance  of  different  sample  times  is  then  easily  determined  from  the  root- 
locus  plot . 


To  establish  transient  response  based  on  system  poles  in  the  z-plane,  the 
relationships 

s = a + jw 

(25) 

are  employed.  Equation  (25)  shows  that 

121  * 

$ • iMI  (26) 

Therefore,  for  a given  At,  (1)  lines  of  constant  a map  into  circles  and  (2)  lines 
of  constant  u map  into  rays  originating  at  the  origin.  To  establish  lines  of 
constant  damping  (£)  in  the  z-plane  the  relationships 

<r  * $wn 


15  1 


(27) 
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along  with  Eq.  (26)  are  used.  Lines  of  constant  damping  (for  £ «■  0 to  1.0)  as 
well  as  lines  of  constant  cAt  and  wAt  are  shown  in  Fig.  3. 

The  system  will  be  stable  for  all  values  of  oAt  for  which  the  root  locus  is 
inside  the  unit  circle  (oAt  < 0).  For  transient  response.  Fig.  3 shows  that: 

(1)  for  a pole  on  the  positive  real  axis  (to  - 0)  and  inside  the  unit  cir- 
cle the  time  response  is  exponentially  decaying.  The  rate  of  decay 
increases  as  the  pole  approaches  the  origin.  There  is  no  decay  if  the 
pole  is  on  the  unit  circle. 

(2)  for  complex  poles  in  the  right -half  plane  and  inside  the  unit  circle 
the  time  response  is  an  exponentially  decaying  sinusoidal  oscillation. 
The  frequency  of  oscillation  increases  as  the  poles  move  toward  the 
imaginary  axis.  The  rate  of  decay  is  faster  near  the  origin.  Stable 
oscillations  result  if  the  poles  are  on  the  unit  circle. 

(3)  for  a pole  on  the  negative  real  axis  (u>At  = n)  as  well  as  for  complex 
poles  in  the  left-half  plane  (and  inside  the  unit  circle)  the  time 
response  is  an  exponentially  decaying  sinusoidal  oscillation.  The  fre- 
quency of  .oscillation  is  higher  than  the  right-half  plane  frequencies. 
Decay  is  faster  near  the  origin. 

Transient  response  for  a discrete-time  system  as  a function  of  z-plane  pole  loca- 
tion is  summarized  in  Fig.  4. 

The  z-transform  for  the  discrete  controller  to  be  implemented  on 
microprocessors  (Eq.  (23))  is  given  by 

u(/)*  K(/)ll£)  (28) 


*<?>*  6»<y-  +D>'’  Ho*. 

Note  that  the  transfer  matrix  K (J)  is  a function  of  sample  time  (see  Eqs.  (23) 
and  (24)).  To  relate  sample  time  (At)  and  system  performance  a root -locus  as  a 
function  of  sample  time  is  plotted.  That  is,  the  poles  of  K (J)  as  a function  of 
At  are  plotted  in  the  s-plane.  Accuracy  of  the  gains  in  the  transfer  matrix  K 
are  limited  by  the  computer  word  length  established  by  applying  the  performance 
index  discussed  previously.  From  the  root-locus  plot  the  maximum  sample  time  for 
which  system  performance  is  adequate  may  be  determined. 
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By  examining  z-plane  poles  as  a function  of  At  the  maximum  sample  time  for 
stable  controller  operation  (all  poles  within  the  unit  circle)  may  be  established. 
In  addition,  controller  transient  response  as  a function  of  sample  time  may  be 
investigated  (see  Fig.  4).  The  maximum  sample  time  for  which  all  controller  real 
poles  are  in  the  right-half  plane  (l.e.,  high  oscillations  associated  with  real 
poles  in  the  left-half  plane  are  not  present  in  the  controller  time  response)  may 
be  determined.  Based  on  this  analysis  the  control  designer  may  determine  the 
maximum  sample  time  required  for  adequate  controller  response.  This  sample  time 
must  be  smaller  than  the  closed-loop  system  sample  time  based  upon  the  Nyquist 
criteria.  Also,  the  sample  time  selected  must  be  long  enough  to  complete  all 
arithmetic  and  I/O  computations. 

Computation  Time 

The  computation  time  per  sampling  interval  includes  (1)  arithmetic  compu- 
tation time  and  (2)  I/O  computation  time.  Arithmetic  computation  time  depends 
upon  (1)  the  number  of  arithmetic  operations  per  sampling  interval  and  (2)  the 
microprocessor  time  to  execute  each  operation.  The  I/O  computation  time  depends 
upon  (1)  the  number  of  inputs  and  outputs  and  (2)  the  speed  of  the  interfaces. 

In  this  study,  it  is  assumed  that  the  interfaces  are  A/D  and  D/A  converters. 

To  reduce  the  number  of  computational  operations,  Eq.  (23)  rather  than 
Eqs.  (13)  and  (14)  are  coded  on  the  microprocessor.  Note  that  the  matrix  T 
should  be  selected  so  that  is  the  desired  structure. 


The  software  to  implement  Eq.  (23)  on  a microprocessor  consists  of  matrix/ 
vector  multiplications  and  one  vector  addition.  For  computational  speed  the 
matrix/vector  multiplications  (p  (q  x 1)  - M(q  x r)  v (r  x 1))  are  accomplished 
by  (1)  multiplying  all  q elements  in  the  jth  column  of  the  matrix  M by  Vj  (i.e., 
Pj,  - M..  v,  for  i - 1 to  q)  for  J ■ 1 to  r and  (2)  adding  the  r columns  of  the 
itn  row  of  P to  obtain  p^  for  i ■ 1 to  q.  This  procedure  for  matrix/vector 
multiplication  is  more  efficient  than  directly  coding  p±  « My  v^  because  Vj 

only  has  to  be  addressed  once  rather  than  r times  for  the  matrix/vector  multipli- 
cation. 


A block  diagram  of  the  matrix/vector  multiplication  code  is  shown  in  Fig.  5. 
Figure  5 indicates  that,  in  addition  to  the  time  for  multiplying  or  adding  the 
appropriate  numbers,  there  is  logic  time  associated  with  each  multiplication  and 
addition.  The  multiplication,  multiplication  logic,  addition,  and  addition  logic 
times  are  denoted  by  M,  A,  Lal,  and  L#2,  respectively.  These  times  are 

a function  of  microprocessor. 
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The  arithmetic  computations  per  sampling  interval  for  the  three  structures  — 
standard,  Jordan  canonical,  and  companion  forms  — are  compared  in  Table  I.  The 
structure  employed  within  the  filter  affects  only  the  filter  equation.  Table  1 
indicates  that  for  n > 1 the  Jordan  canonical  form  requires  the  minimum  arithmetic 
computations  per  sampling  interval.  Also,  the  effect  of  structure  on  computa- 
tional requirements  is  most  significant  when  the  number  of  states  is  significantly 
greater  than  the  number  of  outputs  and  controls  since  differences  in  the  filter 
equation  computational  requirements  are  a function  only  of  the  number  of  states  (n) 

The  computation  time  as  a function  of  addition,  multiplication,  logic,  and 
A/D  and  D/A  conversion  times  as  well  as  state,  input,  and  output  orders  is  shown 
in  Table  I.  The  arithmetic  and  I/O  operation  times  vary  with  microprocessor  and 
A/D  and  D/A  converters,  respectively.  The  characteristics  of  a representative 
set  of  microprocessors  and  A/D  and  D/A  converters  are  shown  in  the  Appendix 
(Tables  A-I  and  A-II,  respectively).  The  microprocessor  characteristics  include 
word  length,  internal  registers,  indexed  addressing  capabilities,  and  multiply 
instruction  capability.  The  A/D  and  D/A  converter  characteristics  include  word 
length,  conversion  time,  and  technology. 

To  compute  the  arithmetic  operation  time, code  must  be  written  to  executed 
the  matrix/vector  multiplication  (see  Fig.  5).  This  code  depends  upon  the  micro- 
processor. Two  microprocessors  were  selected  to  (1)  demonstrate  the  implementa- 
tion of  LQG  control  logic  on  microprocessors  and  (2)  define  the  computation  time 
as  a function  of  system  order  only.  An  Intel  8080  and  a Digital  Equipment  Corpora- 
tion (DEC)  LSI  11/2  — two  very  different  microprocessors  — were  chosen.*  Both  pro 
cessors  are  commercially  available  and  extensively  used.  The  Intel  8080  is  an 
8-bit  microprocessor  with  a limited  instruction  set  and  no  indexed  addressing; 
whereas,  the  DEC  LSI  11/2  is  a 16-bit  microprocessor  with  a powerful  instruction 
set  and  indexed  addressing.  Table  II  lists  representative  microprocessor 
instructions  for  the  Intel  8080  and  DEC  LSI  11/2.  Note  that  hardware  multiply/ 
divide  instructions  are  available  for  the  DEC  LSI  11/2.  * 

The  codes  to  implement  LQG  control  logic  on  an  Intel  8080  and  DEC  LSI  11/2 
have  been  included  in  the  Appendix,  Figs.  Al  and  A2,  respectively.  These  codes 
are  for  the  controller  with  the  standard  structure  in  the  filter.  They  would 
have  to  be  modified  to  take  advantage  of  reductions  in  computational  requirements 
due  to  a Jordan  canonical  or  companion  structure  in  the  filter.  These  codes  use 
a software  multiply  subroutine.  If  a hardware  multiply  were  available,  the  hard- 
ware multiply  would  replace  the  call  to  the  software  multiply  subroutine.  Charac- 
teristics of  hardware  multipliers  are  summarized  in  the  Appendix  (Table  A-III). 
Hardware  multiplier  characteristics  include  word  length,  multiply  time,  and 
technology. 
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Table  III  shows  times  associated  with  executing  the  matrix/vector  multiplica- 
tion on  an  Intel  8080  and  DEC  LSI  11/2.  The  add  and  multiply  logic  times  as  well 
as  the  add  time  are  significantly  less  for  the  DEC  LSI  11/2  than  for  the  Intel  8080. 
This  decrease  in  execution  time  with  an  Increase  in  word  length  is  due  primarily 
to  the  more  powerful  instruction  set  and  addressing  modes  of  the  DEC  LSI  11/2. 

On  the  other  hand,  the  software  multiply  time  for  the  DEC  LSI  11/2  is  4.25  times 
the  execution  time  for  the  Intel  8080.  This  increase  is  due  to  the  fact  that  the 
word  length  of  the  DEC  LSI  11/2  (16  bits)  is  double  the  word  length  of  the  Intel 
8080  (8  bits).  For  comparison,  the  execution  time  for  a double  precision  (16  bit) 
software  multiply  on  the  Intel  8080  is  shown  in  Table  III.  This  time  is  greater 
than  the  software  multiply  time  on  the  DEC  LSI  11/2.  Hardware  and  software 
multiply  instruction  times  for  the  DEC  LSI  11/2  are  also  compared  in  Table  III. 

This  comparison  indicates  that  the  hardware  multiply  instruction  time  on  the 
DEC  LSI  11/2  is  approximately  3.5  percent  of  the  software  multiply  time. 

Intel  8080  and  DEC  LSI  11/2  arithmetic  computation  times  with  software 
multiply  subroutines  are  compared  in  Table  IV.  Table  IV  shows  that  the  computa- 
tion time  per  sampling  interval  for  the  DEC  LSI  11/2  is  significantly  greater  than 
the  computation  time  per  sampling  interval  on  the  Intel  8080  due  to  the  greater 
DEC  LSI  11/2  multiply  time.  Table  IV  also  compares  DEC  LSI  11/2  software 
multiply  arithmetic  computation  time  with  DEC  LSI  11/2  hardware  multiply  arithmetic 
computation  time.  This  comparison  indicates  that  the  arithmetic  computation 
time  per  sampling  interval  on  the  DEC  LSI  11/2  can  be  decreased  by  approximately 
90  percent  by  using  a hardware  multiply.  Table  IV  also  shows  that  computation 
time  is  more  influenced  by  the  number  of  states  (n)  than  by  the  number  of  inputs 
(m)  or  outputs  (£). 

Computation  time  per  sampling  interval  as  a function  of  microprocessor,  system 
structure  within  the  controller,  and  system  order  is  shown  in  Fig.  6.  It  is 
assumed  that  the  system  has  the  same  number  of  states,  Inputs,  and  outputs. 

Figure  6 shows  that  for  a given  sampling  interval  the  number  of  states  that  can 
be  in  the  controller  varies  with  microprocessor.  For  example,  for  a sampling 
interval  of  25  msec  the  maximum  number  of  states  in  a standard  structure  controller 
is:  2 states  for  the  DEC  LSI  11/2  with  software  multiply,  4 states  for  the  Intel 
8080  with  software  multiply,  and  8 states  for  the  DEC  LSI  11/2  with  hardware 
multiply.  Figure  6 also  shows  that  the  computation  time  with  the  Jordan  canonical 
structure  within  the  controller  is  significantly  less  than  the  computation  time 
with  the  standard  structure  within  the  controller. 

The  arithmetic  and  I/O  computation  requirements  per  sampling  Interval  (Table 
I)  establish  the  minimum  sampling  interval.  If  this  minimum  sampling  interval 
is  less  than  the  sample  time  required  to  achieve  adequate  system  performance,  then 
the  mlcroprocessor/interf ace  system  may  be  used  to  implement  the  LQG  controller. 
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Memory  Requirements 

Memory  requirements  depend  upon  (1)  the  system  model  and  (2)  the  computer 
code  including  temporary  storage  to  implement  the  LQG  control  algorithm.  System 
model  memory  requirements  are  a function  of  model  structure  as  well  as  system 
state,  input,  and  output  orders.  System  model  requirements  will  not  vary  with 
microprocessor.  On  the  other  hand,  the  computer  code  and  temporary  storage 
requirements  will  vary  with  microprocessor  as  well  as  system  model. 

Memory  requirements  as  a function  of  system  order  and  structure  are  shown 
in  Table  V for  the  Intel  8080  and  DEC  ISI  11/2.  Table  V indicates  that  the  system 
and  temporary  storage  memory  requirements  are  primarily  dependent  on  the  number 
of  states  in  the  system.  The  standard  structure  requirements  increase  with  the 
square  of  the  number  of  states;  whereas,  the  Jordan  canonical  and  companion  struc- 
tures increase  linearly,  with  the  number  of  states.  Table  V shows  that  memory 
requirements  for  the  Jordan  canonical  and  companion  structures  are  the  same.  The 
Intel  8080  computer  code  memory  requirements  (341  words)  are  approximately  triple 
the  DEC  LSI  11/2  computer  code  memory  requirements  (113  words).  This  difference 
in  code  requirements  is  due  to  the  fact  that  the  DEC  LSI  11/2  instruction  set  is 
more  powerful  than  the  Intel  8080  instruction  set. 

Figure  7 shows  PROM  and  RAM  memory  requirements  as  a function  of  model  struc- 
ture within  the  controller  for  the  Intel  8080  microprocessor.  It  is  assumed  in 
Fig.  7 that  the  number  of  system  states,  inputs,  and  outputs  are  equal.  Figure  7 
indicates  that  memory  requirements  for  the  LQG  controller  are  not  severe.  For 
example,  if  the  system  has  10  states,  10  inputs,  and  10  outputs  then  the  PROM  and 
RAM  memory  requirements  are  0.641  Vi  and  0.365  bytes,  respectively.  On  the  other 
hand,  the  arithmetic  computation  time  per  sample  interval  for  this  10  state  con- 
troller implemented  on  the  PEC  LSI  11/2  with  hardware  multiply  (the  fastest  micro- 
processor considered)  would  be  33,5  msec.  Therefore,  computation  requirements 
are  a more  critical  consideration  than  memory  requirements. 
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APPLICATION  AND  EVALUATION  OF 
MICROPROCESSOR  REQUIREMENT  PROCEDURES 

The  application  and  evaluation  of  microprocessor  requirement  procedures  are 
discussed  in  this  section.  To  evaluate  and  illustrate  the  developed  procedures 
they  were  applied  to  (1)  a second-order  plant  and  (2)  a fifth-order  F100  turbofan 
engine  model  linearized  at  sea  level  static  military  operation.  To  verify  the 
results  a continuous  system/microprocessor  controller  was  simulated  on  the  UNIVAC 
1110  at  UTRC. 

The  closed-loop  simulation  for  verifying  the  procedures  is  described  first. 
Use  of  these  procedures  to  define  microprocessor  requirements  for  digital  control 
of  a second-order  system  are  presented  next.  In  the  final  section  results  of 
applying  the  procedures  to  establish  microprocessor  requirements  of  a linear 
quadratic  Gaussian  (LQG)  controller  for  the  fifth-order  F100  turbofan  model  are 
set  forth. 


Closed-Loop  Simulation  for  Verifying  Procedures 

The  closed-loop  simulation  to  be  used  in  verifying  the  procedures  for 
establishing  microprocessor  requirements  to  implement  l/)G  control  logic  is  shown 
in  Fig.  8.  The  simulation  consists  of  (1)  the  open-loop  system  dynamics  (Eq.  (1)) 
and  (2)  discrete  controller  dynamics  (Eq.  (23)).  Th^  open-loop  system  was 
simulated  using  floating-point  arithmetic  and  a very  small  integration  step 
size  (0.001  seconds)  so  that  the  plant  appears  as  a continuous  system.  The 
discrete  control  and  estimation  logic  were  coded  to  simulate  fixed-point  arithmetic 
operation.  Coding  was  developed  in  the  controller  so  that  the  word  length  and 
sample  time  can  be  varied.  The  word  length  can  be  varied  from  1 to  36  bits. 

The  resulting  closed-loop  system  appears  as  a continuous  plant/discrete  controller. 

The  analysis  for  establishing  the  microprocessor  requirement  procedures  is 
based  on  linear  system  theory-  however,  when  the  IQG  control  algorithm  is  coded 
on  a digital  computer  the  algorithm  is  nonlinear  due  to  finite  word  lengths. 

That  is,  the  matrix  multiply 


p,  = m„  v,  + m,2  v2 


(29) 
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is  calculated  in  the  computer  as 


(P,)c  = 


[^llCV|c)c+  (ml2C  V2c)c]c. 


(30) 


However,  in  the  analysis  it  is  assumed  that 


(p,)c  = (m))c  v,c  + m^c  V2c)c. 


(31) 


This  assumption  in  the  analysis  (Eq.  (31))  must  be  taken  into  account  when  inter- 
preting results  for  very  small  word  lengths  (i.e.,  the  filter  is  very  nonlinear  - 
(Eq.  30)). 

Second-Order  System 

Normalized  plant,  control,  and  filter  dynamics  (Eqs.  (1)  and  (2))  are  shown 
in  Table  VI.  In  control  design  the  deterministic  control  weighting  matrices 
were  selected  as  Identity  matrices. 

The  performance  index  for  the  second-order  plant  as  a function  of  word 
length  and  model  structure  within  the  controller  is  shown  in  Fig.  9.  For 
notatlonal  convenience  the  performance  index  is  normalized  so  that  the  perfor- 
mance index  scale  ranges  from  0 to  1.0.  Figure  9 indicates  that  closed-loop 
system  performance  Improves  significantly  (J  decreases  rapidly)  as  the  controller 
word  length  increases  from  3 to  8 bits.  If  less  than  3 bits  are  employed  in 
the  controller  the  closed-loop  systmn  is  unstable.  As  the  controller  word 
length  is  Increased  from  8 to  16  bits  small  change  in  closed-loop  system  performance 
occurs.  In  addition.  Fig.  9 indicates  that  the  model  structure  employed  within 
the  filter  has  little  effect  on  the  closed-loop  system  sensitivity  to  controller 
word  length  for  word  lengths  greater  than  or  equal  to  seven  bits.  However,  for 
word  lengths  less  than  seven  bits  and  for  this  closed-loop  system  the  Jordan 
canonical  form  is  slightly  less  sensitive  than  the  companion  and  standard  forms 
to  controller  errors  due  to  finite  microprocessor  word  length. 

Output  response  of  the  second-order  model  with  different  word-length  con- 
trollers is  compared  in  Fig.  10.  The  standard  structure  was  employed  In  the 
controller.  Output  response  with  different  structures  and  6 bit  word  length  in 
the  controller  are  shown  in  Fig.  11.  The  controller  sample  time  in  Figs.  10  and 
11  is  0.1  seconds.  The  response  in  Figs.  10  and  11  results  from  initial  condi- 
tions of  0.5  on  the  normalized  states.  Transient  response  of  Figs.  10  and  11 
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verify  results  obtained  from  analytically  computing  the  performance  index  J. 

In  addition,  second-order  system  time  responses  with  different  sample  times  in 
the  controller  were  generated.  Different  controller  sample  times  did  not  change 
the  relative  effect  of  different  controller  word  lengths  on  system  performance. 

Note  that  the  steady-state  bias  errors  shown  in  Figs.  10  and  11  for  controller 
word  lengths  less  than  or  equal  to  8 bits  do  not  affect  the  performance  index 
(Gq.  (17)).  (If  they  did,  the  value  of  the  performance  index  would  approach 
infinity.)  This  result  is  due  to  the  fact  that  the  performance  index  is  com- 
puted from  Eq.  (18);  whereas,  response  of  the  continuous  system/discrete 
controller  simulation  is  computed  from  Eqs.  (1),  (2),  and  (23). 

To  determine  the  required  sample  time  for  the  second-order  8 bit  controller 
a root  locus  of  the  controller  poles  as  a function  of  sample  time  is  plotted. 

Recall  that  the  requirement  for  8-bit  accuracy  is  established  using  the  perfor- 
mance index.  The  root  locus  plot  is  shown  in  Fig.  12.  Poles  at  selected 
sample  times  (0.1,  0.7  and  1.3  seconds)  are  displayed.  Figure  12  indicates 
that  for  sample  times  less  than  or  equal  to  0.7  seconds  the  poles  are  in  the 
right-half  plane  and  are  heavily  damped  (see  Fig.  3).  Therefore,  for  sample 
times  less  than  or  equal  to  0.7  seconds  the  time  response  of  the  controller 
should  not  be  oscillatory.  In  addition,  response  decay  time  will  decrease  as 
the  sample  time  increases  from  0.1  to  0.7  seconds.  Figure  12  shows  that  one  of 
the  poles  is  close  to  the  unit  circle  in  the  left-half  plane  when  the  sample 
time  is  equal  to  1.3  seconds.  Therefore,  for  a sample  time  of  1.3  seconds  the 
controller  response  will  be  stable  but  oscillatory.  The  maximum  frequency  of 
the  closed-loop  system  is  0.54  hertz.  To  satisfy  the  Nyquist  rate  the  sample 
time  must  be  less  than  0.93  seconds.  Therefore,  the  sample  time  of  0.7  seconds 
satisfies  the  Nyquist  rate.  The  maximum  sample  time  for  the  second-order 
controller  is  0.7  seconds.  With  this  sample  rate  the  system  response  will  not 
be  oscillatory. 

Second-order  output  response  as  a function  of  sample  time  is  shown  in  Fig. 

13.  The  standard  structure  with  8 bit  accuracy  was  employed  within  the  controller. 
The  output  response  results  from  an  initial  condition  of  0.5  on  both  normalized 
states.  Figure  13  verifies  results  obtained  from  the  z -plane  plot  analysis. 

Computation  times  required  per  sampling  Interval  for  the  second-order  system 
using  software  multiplication  are  listed  in  Table  VII.  This  table  Indicates 
that  the  minimum  computation  time  is  2.73  msec  for  the  Jordan  canonical  structure 
coded  on  the  Intel  8080.  The  maximum  computation  time  is  10.12  msec  for  the 
standard  structure  coded  on  the  DEC  LSI  11/2.  Both  times  are  well  within  the 
maximum  required  sample  time  of  700  msec. 
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Memory  requirements  for  the  second-order  controller  are  shown  in  Table  VIII. 
Table  VIII  shows  that  model  structure  within  the  second-order  controller  has 
little  effect  on  memory  requirements  due  to  the  low  system  order. 

In  summary,  results  Indicate  that  the  Intel  8080  with  software  multiply  can 
be  used  to  Implement  the  LQG  control  law  for  this  second-order  system.  The 
requirements  for  the  second -order  controller  are  (1)  8-bit  word  length  accuracy, 
(2)  a minimum  sample  time  (based  on  computations)  of  2.73  msec  and  a maximum 
sample  time  (based  on  controller  poles)  of  700  msec,  and  (3)  347  words  of  PROM 
and  21  words  of  RAM. 


Fifth-Order  Gas  Turbine  Engine  Model 

Closed-loop  engine  dynamics  are  shown  in  Table  IX.  The  plant  is  a fifth-order 
F100  turbofan  engine  model  linearized  at  sea  level  static  military  operation. 

The  deterministic  control  design  matrices  were  identity  matrices.  The  process 
and  measurement  noise  statistics  are  shown  in  Table  IX. 

The  normalized  performance  index  for  the  fifth-order  engine  model  as  a 
function  of  word  length  is  displayed  in  Fig.  14.  The  standard  model  structure 
was  employed  within  the  controller  for  the  results  illustrated  in  Fig.  14.  Trans- 
forming the  closed-loop  standard  engine  model  to  the  companion  or  Jordan  canonical 
structures  resulted  in  a numerically  ill-conditioned  closed-loop  matrix.  That 
is,  the  magnitudes  of  the  elements  in  the  closed-loop  matrix  varied  over  a wide 
range.  For  example,  the  magnitude  of  the  companion  form  closed-loop  matrix  elements 
varied  from  0.076  to  5.357  x 10®.  As  a result  of  this  numerical  variation  the  closed- 
loop  system  dynamics  were  unstable  for  word  lengths  less  than  or  equal  to  16  bits. 
Figure  14  shows  that  closed-loop  engine  performance  with  the  standard  model  structure 
within  the  controller  Improves  as  the  word  length  increases  from  7 to  12  bits. 
Increasing  the  word  length  from  12  to  16  bitB  does  not  significantly  affect  closed- 
loop  performance. 

Afterburner  pressure  response  of  the  engine  model  — representative  of 
closed-loop  engine  response  — with  different  word-length  controllers  is  compared 
in  Fig.  15.  The  responses  in  Fig.  15  result  from  Initial  conditions  of  0.1  on 
the  normalized  engine  states.  Figure  15  verifies  results  obtained  from 
analytically  computing  the  performance  index  J.  However,  the  technique  to 
analytically  and  in  closed  form  evaluate  the  performance  degradation  due  to  a 
finite  word  length  in  the  controller  (Eqs.  (21)  and  (22))  eliminates  the  need 
to  perform  these  numerical  integrations. 

To  establish  the  maximum  sample  time  the  poles  of  the  12-blt  discrete 
controller  as  a function  of  sample  time  are  plotted  in  the  z-plane.  Poles 
at  selected  ssmple  times  (0.010,  0.025  and  0.046  seconds)  are  displayed  in 
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Fig.  16.  Figure  16  Indicates  that  for  sample  times  less  than  or  equal  to  0.025 
seconds  the  poles  of  the  discrete  controller  are  in  the  right-half  plane.  Response 
for  sample  times  less  than  or  equal  to  0.025  seconds  will  be  non-oscillatory. 

Figure  16  shows  that  for  a sample  time  of  0.046  seconds  controller  response  will 
be  oscillatory.  For  sample  times  greater  than  0.047  seconds  there  is  a pole 
outside  the  unit  circle  (that  is,  the  controller  is  unstable).  The  maximum 
frequency  of  the  closed-loop  system  is  5.4  hertz.  Therefore,  the  Nyquist  sample 
time  is  0.09  seconds.  The  maximum  sample  time  of  0.025  seconds  is  less  than  the 
Nyquist  sample  time.  Therefore,  the  arithmetic  and  1/0  computation  time  per 
sampling  interval  must  be  less  than  or  equal  to  25  msec  to  achieve  satisfactory 
performance.  Figure  17  verifies  results  of  the  z-plane  plot  analysis. 

Computation  times  required  per  sampling  Interval  for  the  fifth-order  engine 
model  are  listed  in  Table  X.  The  times  shown  are  for  the  standard  structure 
within  the  controller.  Recall  that  controllers  with  the  companion  or  Jordan 
canonical  structure  were  not  stable  for  word  lengths  of  16  bits  or  less.  Also, 
only  computation  times  for  the  DEC  LSI  11/2  microprocessor  are  shown  since  a 
word  length  of  12  bits  is  the  accuracy  required  to  Implement  the  controller. 

Table  X indicates  that  a sampling  interval  time  of  87.48  msec  is  required  to 
complete  the  control  arithmetic  and  I/O  computations  if  a software  multiply 
subroutine  is  used.  Since  the  maximum  sample  time  from  the  z-plane  plot  analysis 
is  25  msec,  the  DEC  LSI  11/2  with  software  multiply  is  not  fast  enough  for 
implementing  the  controller.  On  the  other  hand,  the  total  computation  time 
for  the  DEC  LSI  11/2  with  a hardware  multiply  is  9.68  msec.  This  time  is  less 
than  the  required  maximum  sample  time  of  25  msec.  Therefore,  the  DEC  LSI  11/2 
with  the  hardware  multiply  package  can  be  used  to  Implement  this  LQG  F100  engine 
controller. 

Memory  requirements  (see  Table  XI)  for  the  linearized  F100  engine  controller 
are  less  than  200  words  of  RAM  and  200  words  of  PROM.  These  memory  requirements 
are  very  small. 

In  summary,  the  DEC  LSI  11/2  with  a hardware  multiply  can  be  used  to 
Implement  the  LQG  F100  engine  controller.  The  requirements  for  the  linearized 
F100  engine  model  LQG  controller  are  (1)  12-bit  word  length  accuracy,  (2)  a 
minimum  sample  time  (based  on  computations)  of  9.68  msec  and  a maximum  sample  time 
(based  on  controller  poles)  of  25.0  msec,  and  (3)  less  than  200  words  of  RAM 
and  200  words  of  PROM. 
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CURRENT  UTRC  PROGRAM 


The  current  UTRC  program  sponsored  by  APOSR  Is  a logical  extension  of  the 
research  reported  here.  The  program  Is  directed  toward  Cl)  validating  results 
presented  In  this  report  by  demonstrating  microprocessor  implementation  of  LQG 
control  logic  and  (2)  extending  the  analysis  of  microprocessor  requirements  to 
nonlinear  systems. 

To  demonstrate  microprocessor  Implementation  of  LQG  control  logic,  a system 
that  consists  of  an  analog  computer,  a microcomputer,  A/D  and  D/A  converters, 
and  a display  unit  is  being  employed.  Continuous  linear  system  dynamics  are  simulated 
on  the  analog  computer.  Discrete  LQG  control  logic  is  coded  on  the  microcomputer. 
General  purpose  A/D  and  D/A  converters  are  employed  to  interface  the  microcomputer 
controller  with  the  analog  computer  simulation. 

To  extend  the  analysis  to  nonlinear  systems,  scheduled-gain  nonlinear 
estimation  and  control  logic  based  on  LQG  theory  is  being  employed.  Require- 
ments — including  accuracy,  computational  capability,  and  memory  requirements  — 
for  Implementing  nonlinear  estimation  and  control  logic  are  being  established. 

The  procedures  employed  to  establish  microprocessor  requirements  for  linear 
systems  are  being  extended  to  determine  microprocessor  requirements  for  non- 
linear systems. 
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B 
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C 

C 

D 

E 

Et 
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F 

F 


F^ 


G 

G 

<t> 

G? 


LIST  OF  SYMBOLS 

Constant  n x n matrix  used  to  describe  linear  system  dynamics 
Coefficients  of  the  characteristic  polynomial  of  the  matrix  F 
Constant  n x m matrix  used  to  describe  linear  system  dynamics 
Number  of  bits  in  digital  word 

Constant  p x n matrix  used  to  describe  linear  system  dynamics 

Constant  p x 2n  matrix  used  to  describe  closed-loop  linear  system 
dynamics 

Constant  p x m matrix  used  to  describe  linear  system  dynamics 
Constant  I x n matrix  used  to  describe  linear  system  dynamics 
Constant  I x n matrix  used  to  describe  transformed  linear  system  dynamics 
Quantization  * error 

Constant  n x n matrix  used  to  describe  optimal  deterministic  closed- 
loop  system  dynamics 

Constant  2n  x 2n  matrix  used  to  describe  optimal  stochastic  closed- 
loop  system  dynamics 

Constant  n x n matrix  used  to  describe  transformed  optimal  closed-loop 
system  dynamics 

Constant  m x n optimal  deterministic  closed-loop  feedback  gain  matrix 

Constant  m x 2n  optimal  stochastic  closed-loop  feedback  gain  matrix 

Constant  m x n optimal  deterministic  closed-loop  feedback  gain  matrix 
for  microprocessor  implementation 

Constant  m x n optimal  closed-loop  feedback  gain  matrix  for  transformed 
system 


29 





R79-944258-2 


LIST  OF  SYMBOLS  (Continued) 


Constant  n x I Kalman  filter  gain  matrix 

Constant  n x 1 Kalman  filter  gain  matrix  for  microprocessor  Implementation 

Constant  n x f Kalman  filter  gain  matrix  for  transformed  system  state 

model 

1 x n row  vector  with  hi  * 1.0;  hi  - 0,  1 - 2,  ...  n 

Identity  matrix 

General  subscript 

Performance  index 

General  subscript 

Transfer  matrix  for  discrete  controller 
Discrete  time 

Dimension  of  system  measurement  vector  z 
Dimension  of  system  Input  vector  u 
Milliseconds 

Dimension  of  system  state  vector  x 
Nanoseconds 

Unknown  matrices  in  the  bilinear  equations  for  computing  the 
performance  Index  J (1  “ 1,2,3) 

Dimension  of  system  output  vector  y 

Constant  p x p output  weighting  matrix  in  quadratic  performance  index 

Constant  m x m input  weighting  matrix  in  the  quadratic  performance 
index 
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LIST  OF  SYMBOLS  (Continued) 

Read/Write  Memory 
Programmable  Read  Only  Memory 
Laplace  transform  operator 
Seconds 
Time,  sec 

Constant  nonsingular  n x n state  transformation  matrix 

m x 1 system  input  vector 

m x 1 optimal  system  input  vector 

n x 1 transformed  system  state  vector 

n x 1 system  state  vector 

2n  x 1 closed-loop  system  state  vector 

Closed-loop  system  state  vector  at  t ■ 0 

p x 1 system  output  vector 

I x 1 system  measurement  vector 

Discrete  transform  operator 

1x1  sensor  noise  vector 

Eigenvalue  of  closed-loop  system  (i  • 1,  2,  ...  n) 
Microseconds 

m x 1 input  (process)  noise  vector 

Real  part  of  complex  frequency  variable  s 

T-1  FT  At  + 

2: 
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LIST  OF  SYMBOLS  (Continued) 

Constant  n x n closed-loop  system  matrix  for  microprocessor  implementation 
u Imaginary  part  of  complex  frequency  variable  s 

u±  n x 1 eigenvector  associated  with  the  ith  eigenvalue 

( )j  jth  component  of  vector  in  parentheses 

( )ij  Element  in  the  ith  row,  jth  column  of  the  matrix  in  parentheses 

( )*  System  response  or  controller  matrix  associated  with  36-bit  controller 
+ 

( ) System  response  or  controller  matrix  associated  with  b-blt  controller 

(b  < 36  bits) 

( ) Matrix  in  parentheses  implemented  on  digital  controller 

(')  Time  derivative  of  quantity  in  parentheses 

d(  ) Differential  of  quantity  in  parentheses 

A(  ) Finite  increment  of  quantity  in  parentheses 

A 

( ) Estimated  value  of  quantity  in  parentheses 

( )“*  Inverse  of  matrix  in  parentheses 

( ) ’ Transpose  of  quantity  in  parentheses 

tr(  ) Trace  operator-trace  of  square  matrix  is  equal  to  sum  of  all 
diagonal  elements  of  the  matrix 

A Equals  by  definition 
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:e  Fig.  5 for  definition  of  A, Lai,  La2,  M 
D denotes  time  for  A/D  conversion  (jjiaec) 
A denotes  time  for  D/A  conversion  (/isec) 
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TABLE  II 


REPRESENTATIVE  MICROPROCESSOR  INSTRUCTIONS 

R • Register 
(R)  « Contents  of  Register 
RP  ■ Register  Pair 
PC  * Program  Counter 

INTEL  8080 : Accumulator  is  Register  When  No  Register  Specified 
M Indicates  H,  L Register  Pair 
DEC  LSI  11/2:  R6  is  Stack  Pointer 


Instruction  Description 


Add  Data  From  Memory  to  Register 


Call  Subroutine  Identified  by  Operand 


Double  Precision  Add 


Decrement  Register  by  One 


Input  Data  from  Location  Identified  by 
Operand  to  Register 


Transfer  Program  Control  to  Location 
Identified  by  Operand 


Load  lamed  late  16-Bit  Data  into  Register 


Move  Data  from  one  Register  to  Another 


Move  Contents  of  Register  to  Mmsory 


Move  Immediate  Data  to  Register 


Move  Immediate  Data  to  Memory 


Logical  OR  Register  with  Register 


Output  Data  to  Location  Idantlfled  by 
Operand 


Place  Top  of  Stack  in  Register 


Place  Register  Contents  on  Top  of  Stack 


Return  from  Subroutine 


One's  Complement  Register 


Arithmetic  Rotate  Register  Right/Left 


Rotate  Register  Right /Left 


Logical  AND  Register  with  Register 


Mnil-Polirt  Multiply*1* 


Find-Point  Mrlde*1* 


Floctlng-Polnt  Add*1* 


Floating-Point  Multiply* 


Floating-Point  Dirld*' 


Instruction 


Length  (Bytes) 


Intel  8080  | DEC  LSI  11/2  \ Intel  Q08O 


DEC  LSI  11/2 


(1)  Extandad  Ax  It  teat  le  Chip,  KEVU,  Faquir  ad 


COMPARISON  OF  INTEL  8080  ARITHMETIC  COMPUTATION  TIME  WITH 
DEC  LSI  11/2  ARITHMETIC  COMPUTATION  TIME 
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are  shown  In  the  Appendix  (Table  A-II) . 


R79-944258-2 


TABLE  VI 

SECOND-ORDER  MODEL  DYNAMICS 


■ 5"  V.-  . •, 


* Intel  8080:  1 word  ■ 8 bits 

DEC  LSI  11/2:  1 word  ■ 16  bits 
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TABLE  IX 


FIFTH-ORDER  FlOO  ENGINE  MODEL  DYNAMICS 


Engine  Model 

Linearized 

at  Sea-Level  Static  Military  Operation 

States 

Outputs 

Controls 

Fan  turbine  Inlet 

temperature  Airflow 

Jet 

exhaust  area 

Main  burner  pressure 

Fan  stability  margin 

Fan 

inlet  guide  vanes 

Fan  speed 

Compressor  stability  margin 

Compressor  variable  vanes 

Compressor  speed 

Thrust 

Main  burner  fuel  flow 

Afterburner  pressure 

High  Turbine 

inlet  temperature 

Matrix 

Matrix  Elements 

-34.013 

-9.303 

12.037 

-2.398 

-1.254 

4.389 

-38.762 

-4.221 

28.480 

14.729 

A 

-4.755 

2.287 

-0.400 

-1.546 

-2.200 

2.046 

1.062 

-0.729 

-2.150 

-0.624 

4.150 

-8.814 

-0.167 

7.477 

1.099 

0.766 

0.546 

-0.813 

17.095 

0.056 

1.341 

7.737 

8.641 

B 

0.156 

-1.176 

-0.416 

2.034 

-0.136 

-0.024 

-0.555 

-0.378 

-4.729 

0.874 

1.617 

0.223 

-0.042 

0.063 

0.013 

-0.054 

1.404 

1.045 

0.092 

-0.060 

-0.028 

-0.050 

c 

0.386 

0.100 

-0.217 

-0.095 

0.305 

-0.326 

-0.458 

0.584 

-0.538 

-0.183 

-0.564 

0.394 

-0.165 

0.394 

1.044 

0.001 

-0.013 

0.002 

-0.015 

-0.003 

-0.013 

-0.044 

D 

-0.043 

0.278 

0.035 

-0.155 

-0.101 

0.281 

0.137 

-0.041 

0.073 

0.047 

-0.091 

0.050 

1.0 

0 

0 

0 

0 

0 

n 

0 

E 

0 

0 

1.0 

0 

0 

0 

0 

1.0 

0 

0 

0 

■HI 

1.0 

E 
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TABLE  ix  (Continued) 


Matrix 

Matrix  Elements 

1.602 

-1.183 

2.224 

0.148 

5.530 

(l 

0.012 

3.074 

-0.341 

-0.903 

-0.223 

-2.942 

-5.064 

5.544 

-2.222 

8.148 

-4.362 

0.749 

-0.652 

-0.092 

-0.811 

6.110 

1.879 

20.560 

2.456 

0.036 

3.341 

1.046 

12.150 

3.728 

0.018 

H 

0.822 

0.273 

3.493 

2.523 

0.006 

0.053 

0.045 

1.357 

3.070 

-0.002 

0.064 

0.018 

0.282 

-.189 

0.001 

0.01 

0 

0 

0 

0 

0 

0.01 

0 

0 

0 

Covariance  of 

0 

0 

0.01 

0 

0 

£ 

0 

0 

0 

0.01 

0 

0 

0 

0 

0 

0.01 

0.01 

0 

0 

0 

0 

0.01 

0 

0 

Covariance  of 

0 

0 

0.01 

0 

n 

0 

0 

0 

0.01 

0 

0 

0 

0 

0.01 

1 
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TABLE  XI 

MEMORY  REQUIREMENTS  FOR  FIFTH-ORDER  ENGINE  MODEL 

Standard  Structure  Within  Controller 
DEC  LSI  11/2  Microprocessor 


MICROPROCESSOR 

IMPLEMENTATION 


mmm 
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FIG.  5 


BLOCK  DIAGRAM  FOR  MATRIX/VECTOR  MULTIPLICATION 

1 M(q  x r)v(r  x 1) 

P«  ‘ M1jvi 


1 

p(q  X 1)  ■ 
P 

\ 

(1VSET  COLUMN 
COUNTER  (j) 

(2)  SET  UP  v, 

2 

(1)  SE1 
COUN 

(2) SET 

r row 

TER  (i) 

UP  M„ 

CALL  SOFTWARE 
MULTIPLY  SUBROUTINE 
(OR  HARDWARE 
IMULTIPLY  INSTRUCTION)! 


(l)SET  COLUMN 
COUNTER  (j) 

(2)  SET  UP  PH 

9 

ADD 

INSTRUCTION 
(P,  -P,  +P||) 

NO 

—/COLUMNSS 

NO 


H - TIME  TO  EXECUTE  BLOCK  3 
L^,  - TIME  TO  EXECUTE  BLOCKS  2,4,5 
Cn2  * TIME  TO  EXECUTE  BLOCKS  1,6 
X - TIME  TO  EXECUTE  BLOCKS 

- TIME  TO  EXECUTE  BLOCKS  8,10 
t^2  - TIME  TO  EXECUTE  BLOCKS  7,11 


kjl-q7)> 


YES 


7S-03-1M-6 
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FIG.  6 


ARITHMETIC  COMPUTATION  TIME 

n ■ m ■ I 

a)  COMPUTATION  TIME  AS  A FUNCTION  OF  MICROPROCESSOR  (msec); 
STANDARD  STRUCTURE  WITHIN  CONTROLLER 


— • INTEL  8080 

— DEC  LSI  11/2 

DEC  LSI  11/2  WITH  HARDWARE  MULTIPLY 


b)  COMPUTATION  TIME  AS  A FUNCTION  OF  STRUCTURE  WITHIN  CONTROLLER  (nwcl; 
DEC  LSI  11/2  MICROPROCESSOR 


NUMBER  OF  STATES 


79-03- IB* -8 


y **■ , » v » * ,• 
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FIG.  11 
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SECOND-ORDER  OUTPUT  RESPONSE  AS  A FUNCTION  OF  MODEL 
STRUCTURE  WITHIN  CONTROLLER 

RESPONSE  TO  INITIAL  CONDITION:  X'Q  - |0.5. 0.5) 

At -0.1  SEC 


DIGITAL  CONTROLLER  WITH: 

— 36  BITS,  JORDAN  CANONICAL  STRUCTURE 

— 6 BITS.  JORDAN  CANONICAL  STRUCTURE 

— 6 BITS.  COMPANION  STRUCTURE 


7*— 03— 46— 6 
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FIG.  13 
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SECOND-ORDER  OUTPUT  RESPONSE  AS  A FUNCTION  OF  SAMPLE  TIME 

RESPONSE  TO  INITIAL  CONDITION:  Xq  - (0.6,0.51 
STANDARD  STRUCTURE  WITHIN  CONTROLLER 


DIGITAL  CONTROLLER  WITH: 

■ ~ 36  BITS,  At  - 0.1  SEC 

— — 8 BITS,,  At -0.1  SEC 

— — — ——  8 BITS.  At  -0.7  SEC 
8 BITS.  At  - 1.3  SEC 
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APPENDIX 

MICROPROCESSOR  SURVEY  AND  COMPUTER  CODE 


Characteristics  of  (1)  microprocessors,  (2)  A/D  and  D/A  converters,  and 
(3)  hardware  multipliers  are  presented  in  this  Appendix.  The  characteristics 
tabulated  here  were  obtained  from  Electrical  Design  News  (EDN)  1976-1978  as  well 
as  from  TRW  product  sheets.  Microprocessor  characteristics  — including  word 
length,  internal  registers,  indexed  addressing  capabilities,  and  multiply 
instruction  capability  — are  listed  in  Table  A-I.  The  A/D  and  D/A  character- 
istics — Including  word  length,  conversion  time,  and  technology  — are  shown 
in  Table  A-II.  Table  A-III  displays  multiplier  characteristics  including  word 
length,  multiply  time,  and  technology. 

In  addition,  microprocessor  code  required  for  calculating  the  matrix/vector 
multiplication  computation  times  and  code  memory  requirements  are  presented  in 
this  Appendix.  Codes  for  implementing  linear  quadratic  Gaussian  control  logic  on 
an  Intel  8080  microprocessor  and  on  a DEC  LSI  11/2  microprocessor  are  shown  in 
Figs.  A-l  and  A-2,  respectively.  The  number  of  instructions  for  the  Intel  8080 
microprocessor  is  more  than  double  the  number  of  instructions  for  the  DEC  LSI  11/2 
due  primarily  to  the  limited  addressing  capability  of  the  Intel  8080  microprocessor. 
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(1)  Intel  8022  only 


Converter 

Type 

Manufacturer 

Model 

Word  Length 
(bits) 

Conversion 
Time  (usee) 

Technology 

TRW 

TDC1007J 

8 

35x10-3 

Bipolar 

TRW 

TDC1001J 

8 

400xl0~3 

Bipolar 

TRW 

TDC1002J 

8 

1 

Bipolar 

Analog  Devices 

AD75705 

8 

40 

CMOS 

Datel 

ADC-MC88C 

8 

500 

Bipolar 

A/D 

Analog  Devices 

AD7570L 

10 

120 

CMOS 

Datel 

ADC-HX12B 

12 

20 

Hybrid 

Analog  Devices 

AD572BD 

12 

25 

— 

Mlcre  Networks 
National 

ADC80 

12 

25 

Hybrid 

Semiconductor 

ADC1210 

12 

50 

Hybrid 

TRW 

TDC1016J 

8 

35xl0"3 

Bipolar 

Analog  Devices 

AD7523JN 

8 

100xl0"J 

— 

Datel 

National 

DAC-UP88 

8 

2 

Bipolar 

Semiconductor 

DAC0800 

8 

135 

Bipolar 

Datel 

DAC-088 

8 

150 

Bipolar 

D/A 

TRW 

TDC1017J 

10 

50xl0-3 

Bipolar 

Analog  Devices 

AD7541KN 

12 

1 



Datel 

Harris 

DAC-HK12B 

12 

3 

Hybrid 

Semiconductor 

Harris 

Hl-5612 

12 

85 

Bipolar 

Semiconductor 

Hl-562 

12 

200 

Bipolar 

Datel 

DAC-HA12B 

12 

500 

Hybrid 

Analog  Devices 

AD7531 

12 

500 

Hybrid 

t! 

B 

0 
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*v~. — — : — m — — — 1 — r 


REPRESENTATIVE  MULTIPLIERS 


J I 


R79— 944258— 2 


SSS^»«“sh*010®  " 


FIG.  A— 1 


PRELIMINARY  INTEL  «080  SOFTWARE  FOR  LOG  CONTROLLER 


SSnrroprti^"^ONs!B?HKsSOM“*RE"PRELI,'1N*',t 

{UTRC l»78 


22: 

ii\ 

IV, 

71: 

JO: 

III 

J5: 

W, 

It: 

IV: 

*3i 


»*: 

ill 

U\ 

s*» 

IV, 

til 

l?,*PLPt 

til 

M! 

Hi 

Hi 

(l! 


8R6cS8°Sohm 

H7  COU  S02CN 

B2se?Hos?SJ3h 

NSS1  too  S020H 
NV  COU  SO JOM 
it 2 COU  SoIcM 

IlflffWItllGN 
NU  COU  50RSH 

ss  m \ 

pc  cou  i 

N2MJ  COU  0 

«S‘cMu.U. 

•CL  COU  SOM 
MV  1 A * MM 
LII  M.STt 
INK  M 
MV1  B t MS 
NOV  C.M 
1 NX  M 
NOV  C 9 N 

nig  k 

PUSH  PSH 

cN*oLvi.:rT 

L X I M * N 1 
NOV  E,H 
NV I D.RCL 
STAX  0 

28!  RsE 

NOV  N.A 
POP  Pin 
POP  H 

882  8 
JN2  NCLP 
PUSH  PS V 
PUSH  H 
LX  I H,NJ 
POV  A , N 

SUI  N?M1 
NOV  N 9 A 
POP  H 
POP  PSM 
OCR  A _ 

JNZ  NPLP 

Hi  ft’"5 

PUSH  • 

WiHc!W 

CALL  SUN 

r«!MH%ss 

NOV  Cs" 

w.l»BCt 

IBS  8 

NOV  N,C 

(88  (sv 


!«iRM^RNGTrS5s8.5!i.TfRp-c-!Sf?tON  *-° 

i po inter  to  phix  partial  rcsults 

iRiBUI  ?8  MRTm  RESULTi 

{POINTER  TO  PHIX 

{POINTER  TO  VECTORS 

{START  ADDR.  FOR  PHIX  CONSTANTS 


{POINTCR  To  XlA«l/A*I» 

128:  8f  RIaIu&cncnts 

{NO.  OF  CONTROLS 


iwavtsraip  .OF  CONSTANTS 

! lHcs?U7SBIu2ESI?5s¥^TgS 


NT  VECTOR  INTO  ACC. 
CSS  OF  CONSTANTS 


{NEXT  ARGUMENT  ADDR. 

: SI 2E  OF  STATE  VECTOR 
{GET  MULTIPLIER 

t gBThul?iplicano 

{STORE  NN.ADOR.  POINTCR, FLAGS  AND  NS 
• 

| Jfiutf  *!n!8E ACC. I TRUNCATE! 

{ ADDR. OF  STORAGE 
{GET  ADDRESS 

{ EXTEND  DC  REG  FOR  AOOR. 

{STORE  RESULT 

; ADO  DISPLACEMENT 

{ R^STO^F* ri AGS, NN, CONSTANTS  LOL ATION9NS 

8 

{test  loop 

{SAVE  FLA6S.ACC. 

{SET  UP  FOR  NEXT  MCA  SURE  HE  NT 


I TEST  OUTER  LOOP 


{COMPLETE  MATRIX-VECTOR  NULTIPlV  IT  ADDIN6 
{ADDITION  SUBROUTINE 
{ PREPARE  FOR  STORING  RESULT 

iilUMBi'tttsiT 

INCH  AOOR. 

(STORE  NCM  AOOR. 

| REST ORE  REGISTERS 


(CONTINUED) 


R79— 944258-2 


mis  pa®  is  best  quality  mcttCABu 

noa  OQPY  fURBISHBD  TO  DDC  — 


FIG.  A— 1 


PRELIMINARY  INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 

(CONTINUED) 


Tfc  t 

?l! 

79s 
• 0: 

81:N»LP1: 

82: 

83: 

89: NCLP1 t 

II! 

37: 

88: 

89: 

90: 

9i: 

92: 

II! 

95: 

II! 

98: 

99: 

100: 

101: 

18!! 

108: 

181! 


119: 

lil.-S-LPl: 

till 


11:9*17): 

39: 

9 0: 

9HNCLP3: 

32: 

193: 

981 

II! 

97: 

R|: 

99: 

SC: 


NVl  A.NS 
LXI  H.ST2 

In*  h 
NVl  8 t NS 
NOV  C,M 
INX  H 

mw 

POSH  N 
PUSH  PSH 
CALL  MULT 
NOV  A,£ 

LX  1 N»N2 
NOV  C.H 

ni*°6RtL 

nov  a.e 

S8i  8f» 

pop  pSv 

POP  M 

pop  * 

OCR  0 

«iN%kS‘ 

PUSH  H 
LX  1 H.N2 
NOV  A.H 
5U1  N?H? 
NOV  N * A 
POP  H 
POP  PSH 
OCR  A 
JN2  NPLPl 
NVl  A.NS 


PUSH  »S* 

NVl  C • NS 
CALL  SUN 
PUSH  M 

nov  e:Sssl 

NVl  Ot RFL 

MS'.8 

INK  0 
NOV  N,C 
POP  H 
POP  PSH 

JN2  SNLP1 
POP  • 

Hi  b:Ks 
WlVsSS"' 

NVl  A.NS 

Hi  R,STJ 

NVl  8 f NC 

MS  S’" 

KShc«" 

PUSH  M 

me  ae« 

NOV  AfC 

Hi  r:S’ 

Ml  D.RCL 


J Hi  pfc88uC7  C0HPLC7C 
; PREP  ARE  FOP  PH1XIK/K)  PRODUCT 
{POINTER  TO  PHI 

jlifENU0[7H£n.¥tCT0R 

5 GE T NULT 1PL1CAN0 


; TCI  • (Cl  = > « OC I 
.•PREPARE  TO  S TOPI  PfSULT 


SSTORC  RESULT 
; REPEAT  OF  LINES  92-97 


iPHIIIK/M  CONPLETC 
{PREPARE  FOR  VECTOR  ADD 

; POINTER  TO  PPEOICT  SOLUTION 

| PREPARE  FOR  6XIR»1/M*1I  NPV 

I NEXT  ARGUNENT  ADDR. 
sCOLUNN  COUNTER 

! lijR^HlAHiJMoiNTCRtFLAOS  ANO  NC 

t 

i ICI9IC >S>IDC> 

1 RESULT  INTO  ACC. (TRUNCA TC I 

itfMMir" 

I EXTEND  DC  *EG  FOR  AOOR. 


(CONTINUED) 


78-03-89—3 


R79— 944258— 2 


SBZ8  PlflQf  IS  B8ST  QUALITY  PRACIICABH 
irm  oorY  ruraiisHSD  to  doc 


F IG.  A— 1 


PRELIMINARY  INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 

(CONTINUED) 


ill: 

STAX  0 

NOV  A,E 

\\i\ 

ISSt 

Ml  N|A 
POP  Piu 

uv. 

&8P  S 

ISO: 

OCR  0 

I5«: 

loot 

101: 

PUSH  M 

02: 

03: 

tiSi  S:S3 

0*: 

16S: 

SUI  N2N3 
HOV  H.A 

100: 

POP  H 

67: 

ill 

POP  PSM 
OCR  A 

109: 

JN2  NRLPS 

170: 

NVI  A.NC 

171: 

INK  H 

172: 

PUSH  B 

l 7 titlil  Pk  : 

PUSH  PSU 

17*: 

HV1  C > NS 

HI: 

m r 

177: 

LXl  H.NU 

ill! 

SS!  oIrel 

HR 

STAX  D 

I NX  D 

102: 

NOV  H,E 

103: 

POP  H 

10*: 

POP  PSW 

10S: 

OCR  A 

tit: 

m rtM 

100: 

JNP  OONE 

\M\ 

191: 

?f:NULT: 

NOV  A . C 

9*: 

ORA  A 

95: 

901 

*8.Tf 

III 

ORA  A 

JN  HUL  3 

199INUL1: 

CALL  INUL 

J3f,NUL2, 

n«  * 

202: NUL3 : 

OCR  E 

mi 

m **c 

205: 

mm\ 

NOV  E ■ A 

K!toW 

m 

"SI  *•* 

210: 

NOV  E , A 

i i: 

m 

1 ’! 

1 SlNULSl 

58 1* 

i R 

[« 

2 •: 

NOV  C,A 

HSi 

8HS' 

1 i* 

JN  RUV* 

Hi: 

tt«"  I.C 

if;.1 

NOV  E*A 

; STORE  RESULT 
I 


tot  uBi'ftmmj.i 

; RE  ST  ORE  FLAGS, NC.CONSTAN 


NSTANTS  LOCATION, NR 


; TEST  LOOP 

; 

: 

! 


J 

I TEST  OUTER  LOOP 

•PREPARE  FOR  ADDITION 

i 

t 

i . 


f EXTEND  FOR  AOOP. 

i arMi:*1  pcsult 

t STORE  NEW  AOOP. 


GX'K*1/K»1)  r UIKM/AMI 

i 

: 

[check  SIGN  OF  ( C » 

j infix 'sHn  of  te» 

1 €£*  = >— *E» 

I PA6ICIPNA6IE  > = > <0E » 

|>*S  CONP.  iei 

i 

U'l'tlU.  ion 

jl*S  CONP.  TOEI 


2*1  CONPlEWENT  OF  1C  I 


8W 

i i-ciPin 

j 2'S  CONP 


MEWm 


in 


(CONTINUED) 


7R-03-RN-1 


R79— 944268— 2 


IBIS  PAGE  IS  BIST  QUALITY  PRACTICABLI 
Jltoa  OOFY  FURBISHED  TO  DDC 


FIG.  A— 1 


[ 


l 


PRELIMINARY  INTEL  8080  SOFTWARE  FOR  LOG  CONTROLLER 

(CONTINUED) 


227:1MUL: 

IIS; 

1«ULI : 

iU\ 

21*: 

i 

iu; 

219: SUM: 
|*£:SUM0: 

2*2: 

iU\ 

2*  S: VSUM : 
g* j:VSU*0: 

211: 

' 1* 


25  61 

III!00"" 


JMP  MUL1 


SSi 


0;0H 


1*1  H,OH 
“ B.6 
H 


s u 

R6L 
JNC  CMNT 

060  0 
OCR  B 
JM2  IMUL1 
!)G 


»sv 


*R6  6 

58°  R 

OCR  C 

ffl  SU"° 

XR*  6 

58$  S 

600  N 

M.Mo 

:« it  1 

MOV  E,6 

385  $SUN 
OCT 


m 


ICLE6R  ML 

ML  LEFT 
S 6CC  LEFT  THRU  C 6RRT 
; CARRY =71 

; YLS|PR00=PR0D*HULTIPL1C6N0 

I 

j PUT  RESULT  IN  Ot 

1 CLF6R  6CC. 

{ 600  I ML  I TO  6CC 
• 

; k ET URN- -NEXT  TIME  INCREMENT 


I8VKE 


TB— 03— M— * 


R 79 -944 258- 2 


IBIS  PA®  IS  SJBST  QUALITY  FRACT1CA&LA 

IBM  OOPY  PUMISHSD  TO  DDC 


FIG.  A— 2 


PRELIMINARY  DEC  LS1 11/2  SOFTWARE  FOR  LOG  CONTROLLER 


i 


5*1 

III 

II 

9*H* 

I?!*  * 

Wt 

j jtNMULl: 

ft! 

IT* 

Ml 

20* 

Ml"11 

UlltOMl 

25* 

|^jL00P2: 

26* 

29* 

I Os 
Sit 

Hi 

J4*pmi* 

gSjXMATs 

37* 

J8S NMUL2: 
39* 

:?i 

42* 

43: 

49: 

4 S:  CT  2: 

46* 

4 4:  LOOPS : 
49: 

5 jtLOOPR* 


5»* 

US 

5v:su«: 

Hi 

62: 

UU: 

65**4: 

Ml 

66:mnuL3: 

69: 

Tut 

Til 

HI 

?SI 


(OCC  LSI  11/2  COOC  FOR 

(KALMAN  FILTER /CONTROLLER —PRELIMINARY 

{DEVELOPED  BY  R.  S.  FIDENS 

JOIRC I9T9 


88$  5s;r*s?*wo 

NOV  SH.R3 


885  Sjfi?0' 


MOV  IR0I»-IR6I 
MOV  ?8ilV^|ROI 

m SSpi 

m am 

NOV  M,RO 

?88  Si"* 

888  H» 

•Nl  L00P2 

m R^,t,,* 

BMC  L00P1 

885  sew* 

MOV  SPMI.RJ 

88!  Si?!"01 
58!  Jcfittf*' 

MOV  IR61* » (ROI 

m $g$°'tooo«*3» 

• 10  C *2 
OEC  2IRM 
BNE  MMUL2 
MOV  5P?,#l 
MOV  SR’,R2 

885  8:88 

CLK  R* 

AOO  (P1|*,R4 

ore  -*3 

»NL  L00P4 
MOV  C«(|62|» 

OEC  41 
bM  LOOPS 
MC  V 191, »1 
MOV  SPP.R? 

MOV  6, M* 

AOO  I»1I*,I42I* 

Slf  SUM 

MOV  POlNtCK.RO 
MOV  KL.IBri 
MC  V SGjpJ 
MOV  N,*'(ROI 

885  *p£|?-(R6I 

JSR  PC. MULT 
MOV  (»G»,(NCI 
MOV  PROP, ICO j I US  I 


Oft  • 

2(iiPI 

gNE^MMULl 


mirw  i°-mu  or  "-**wr 

! START  ADDR.  OF  M-ARRAY  (KALMAN  GAIN$| 

|?l«T°W»-;Bgr2-»RRAY 

{ PUSH  RO 

j^ALJg  MPY., RESULT  IN  R2 

*o8ne*  0r  HZ  wrooucTS 
i YES, CONTINUE 

!88«mt C0UNU" 

! YES, NEXT  ROM 

lilttl  *.BB8:  85  8TOIS 
jtifi  85  5:558188 

iSCALOR  AOO,  RESULTS  IN  R4 
{DONE? 

8 STORE  RESULTS 
{CONE  HITH  M2  PROOUCT 

iB#TS» 

} ST  ART  ADDR,  OF  PMI  (SYSTEM  MATRIXI 

|!»8t°U5':£5508h.t 

1s8Hlor0mpy., RESULT  IN  R2 
}POP  RO 

BtsotTs 

.•  ns 

| S? ART  AbOR.  OF  PRODUCTS 
{START  ADDR.  OF  RESULTS 

1 1151  8f  5:518188 

; SCALOR  ADD.  PESULT  IN  P4 
sSfn&E  RESULT 


5 DONE  h ITH  PHI  XMA  T 
; PREPARE  F0»  VECTOR  A DC 

{VECTOR  AOO 

; DONE  fcITH  VECTOR  ADD? 

{PLRCAT  LINES  T TO  SI 

{ ?V  ?2  T ^ AUOR^£  0^  3 6**  ( CONTFOL  GAINS) 


(CONTINUED) 


R79— 944258-2 


THIS  PASS  IS  BSST  QUALITY  FKACTICABU 
FROM  OOPY  FURBISHED  TO  DOC 


FIG.  A-2 


PRELIMINARY  DEC  LS1 11/2  SOFTWARE  FOR  LOG  CONTROLLER 

(CONTINUED) 


76: CT Is 

MOV 

SP3.RI 

77: 

MOV 

SUl2 

7b: 

MOV 

T9:LC0PS: 

MOV 

N , R 3 

lit:  L00P6 : 

SbB 

R9 

1 R1 ) ♦ . R9 

A 2: 

DEC 

R 3 

U\ 

W 

kS??*2»- 

AS: 

OEC 

RO 

86: 

BNE 

L00P5 

•TsMULT: 

MOV 

IR1M.R2 

88: 

89: 

B M 1 
MOV 

T«^f..R9 

90: 

BMI 

HUL3 

PC, IMUL 
PC 

93: MUL3: 

DEC 

R9 

99: 

COM 

R9 

95 : mulr  : 

JSR 

PC, IMUL 

9b: 

OEC 

R2 

97: 

98: 

18HmuL2 

99:MUl5l 

OEC 

R2 

100: 

CON 

R 2 

101: 

MOV 

<R3I*,R9 

foil 

§hcP 

MULR 

R9 

139: 

COM 

R9 

105: 

BR  MUL1 

10b: IMUL: 

MOV 

A 0, PROD 

107: 

MOV 

VO.PROOl 

§1: ImuL 1 : 

MOV 

ASL 

m«j° 

10: 

ASL 

PROD 

11: 

AOC 

PR001 

112: 

POL 

R2 

11  3: 

bCC 

chnt 

119: 

AOO 

1 R8 1 » PROD 

US: 

tn 

PP001 

m:CHNT: 

RO 

BNE 

1MUL1 

118: 

119: 

RTS 

END 

PC 

; 


U<K«1/K«1)  COMPLETED? 

CHECK  SIGN  MULTIPLIER 

CHECK  SIGN  MULTIPLICAND 

INTEGER.  UNSIGNED  MPT. .RESULT 

2*S  COMPLIMENT  multiplicand 

INTEGER,  UNSIGNED  MPT. .RESULT 
2*S  COMPLIMENT  result 


POSITIVE 

NRONG  S16N 


{2*S  COMPLIMENT  MULTIPLIER 

; CHECK  SIGN  OF  MULTIPLICAND 

{2*S  COMPLIMENT  MULTIPLICAND 

; END  OF  MULTIPLY  MAIN 

{PREPARE  32  BIT  REGISTER  FOR  PnOOUCT 

•$h¥f?£32  BITS  LEFT 
{ADJUST  CARRY 

{CHECK  FOR  1*5  IN  MULTIPLIER 
{ I IN  MULTIPLIER? 


YES .PRODUCT :PROOUCT 
ADJUST  CARRY 

sr* 

YES 


♦ MULTIPLIER 


79—03— 06— 2 


V v.v  . 

1 r r if  I l il  rtf  iUi  II 


Distribution  List 
Report  R7 9-944258-2 
Contract  F49620-78-C-0017 


Director  of  Mathematical  and  Information 

Sciences 

Bldg.  410 

Bolling  AFB,  Washington,  D.  C.  20332 
Attn:  Major  Charles  L.  Nefzger  (16) 

ASD/XRT 

Wright-Pat terson  AFB,  OH  45433 
Attn:  Lt  Col  W.  Othling 

OASD/PA&E 

Pentagon 

Washington,  DC  20360 
Attn:  Mr.  T.  P.  Christie 

NASA-Flight  Research  Center 
Edwards  AFB,  CA  93523 
Attn:  Mr.  H.  A.  Rediess 

SAMSO/YAD 

Technology-Guidance  & Control 
Los  Angeles  AFS,  CA  90045 

NASA-Ames  Research  Center 
Moffett  Field,  CA  94035 
Attn:  Mr.  Elvood  C.  Stewart 

Naval  Ocean  Systems  Center 
Tactical  Command  Control  Division 
San  Diego,  CA  92152 
Attn:  Dr.  Robert  Kolb,  Code  824 

NASA-Langely  Research  Ctr 

Hampton,  VA  23365 

Attn:  Mr.  Larry  W.  Taylor,  Jr. 

AFIT/ENE 

Wright-Pat terson  AFB,  OH  45433 
Attn:  Capt.  Gary  Reid 

Office  of  Naval  Research 
Technology  Group 
800  North  Quincy  St. 

Arlington,  VA  22217 
Attn:  Mr.  David  Siegel 


AFIT/ENE 

Wright-Patterson  AFB,  OH  45433 
Attn:  Capt  James  Negro 

AFFDL/FGL 

Wright-Patterson  AFB,  OH  45433 
Attn:  Paul  Blatt 

RADC/ISCP 

Griffiss  AFB,  NY  13441 
Attn:  Heywood  Webb 

ARO 

P.0.  Box  12,211 

Research  Triangle  Park,  NC  27709 
Attn:  Dr.  J.  Chandra 

Office  of  Naval  Research 
Mathematics  Group 
800  North  Quincy  Street 
Arlington,  VA  22217 
Attn:  Dr.  Stewart  Brodsky 

AFAPL/TBC 

Wright-Patterson  AFB,  OH  45433 
Attn:  Mr.  Charles  Sklra 

AFATL/DLY 

Eglin  AFB,  FL  32542 
Attn:  Dr.  Jesse  Gonzales 

AFWL/LR0 

Klrtland  AFB,  NM  87117 
Attn:  Lt  Col  Dale  Neal 

David  Taylor  Naval  Ship  R&D  Center 
Code  2730 

Annapolis,  MD  21404 

Attn:  Mr.  Walter  J.  Blumberg 

AFWL/AL 

Klrtland  AFB,  NM  87117 
Attn:  Maj  Kenneth  Herring 


